Revealing myopathy spectrum: integrating transcriptional and clinical features of human skeletal muscles with varying health conditions.

Myopathy refers to a large group of heterogeneous, rare muscle diseases. Bulk RNA-sequencing has been utilized for the diagnosis and research of these diseases for many years. However, the existing valuable sequencing data often lack integration and clinical interpretation. In this study, we integrated bulk RNA-sequencing data from 1221 human skeletal muscles (292 with myopathies, 929 controls) from both databases and our local samples. By applying a method similar to single-cell analysis, we revealed a general spectrum of muscle diseases, ranging from healthy to mild disease, moderate muscle wasting, and severe muscle disease. This spectrum was further partly validated in three specific myopathies (97 muscles) through clinical features including trinucleotide repeat expansion, magnetic resonance imaging fat fraction, pathology, and clinical severity scores. This spectrum helped us identify 234 genuinely healthy muscles as unprecedented controls, providing a new perspective for deciphering the hallmark genes and pathways among different myopathies. The newly identified featured genes of general myopathy, inclusion body myositis, and titinopathy were highly expressed in our local muscles, as validated by quantitative polymerase chain reaction.

The functional role of CST1 and CCL26 in asthma development.

BACKGROUND: Asthma is the most common chronic disease in children with an increasing prevalence. Its development is caused by genetic and environmental factors and allergic sensitization is a known trigger. Dog allergens affect up to 30% of all children and dog dander-sensitized children show increased expression of cystatin-1 (CST1) and eotaxin-3 (CCL26) in nasal epithelium. The aim of our study was to investigate the functional mechanism of CST1 and CCL26 in the alveolar basal epithelial cell line A549. METHODS: A549 cells were transfected with individual overexpression vectors for CST1 and CCL26 and RNA sequencing was performed to examine the transcriptomics. edgeR was used to identify differentially expressed genes (= DEG, |log2 FC | ≥ 2, FDR < 0.01). The protein expression levels of A549 cells overexpressing CST1 and CCL26 were analyzed using the Target 96 inflammation panel from OLINK (antibody-mediated proximity extension-based assay; OLINK Proteomics). Differentially expressed proteins were considered with a |log2 FC| ≥ 1, p < .05. RESULTS: The overexpression of CST1 resulted in a total of 27 DEG (1 upregulated and 26 downregulated) and the overexpression of CCL26 in a total of 137 DEG (0 upregulated and 137 downregulated). The gene ontology enrichment analysis showed a significant downregulation of type I and III interferon signaling pathway genes as well as interferon-stimulated genes. At the protein level, overexpression of CST1 induced a significantly increased expression of CCL3, whereas CCL26 overexpression led to increased expression of HGF, and a decrease of CXCL11, CCL20, CCL3 and CXCL10. CONCLUSION: Our results indicate that an overexpression of CST1 and CCL26 cause a downregulation of interferon related genes and inflammatory proteins. It might cause a higher disease susceptibility, mainly for allergic asthma, as CCL26 is an agonist for CCR-3-carrying cells, such as eosinophils and Th2 lymphocytes, mostly active in allergic asthma.

Effects of exercise on whole-blood transcriptome profile in children with overweight/obesity.

BACKGROUND: The current knowledge about the molecular mechanisms underlying the health benefits of exercise is still limited, especially in childhood. We set out to investigate the effects of a 20-week exercise intervention on whole-blood transcriptome profile (RNA-seq) in children with overweight/obesity. METHODS: Twenty-four children (10.21 ± 1.33 years, 46% girls) with overweight/obesity, were randomized to either a 20-week exercise program (intervention group; n = 10), or to a no-exercise control group (n = 14). Whole-blood transcriptome profile was analyzed using RNA-seq by STRT technique with GlobinLock technology. RESULTS: Following the 20-week exercise intervention program, 161 genes were differentially expressed between the exercise and the control groups among boys, and 121 genes among girls (p-value <0.05), while after multiple correction, no significant difference between exercise and control groups persisted in gene expression profiles (FDR >0.05). Genes enriched in GO processes and molecular pathways showed different immune response in boys (antigen processing and presentation, infections, and T cell receptor complex) and in girls (Fc epsilon RI signaling pathway) (FDR <0.05). CONCLUSION: These results suggest that 20-week exercise intervention program alters the molecular pathways involved in immune processes in children with overweight/obesity.

Optimized single-cell RNA sequencing protocol to study early genome activation in mammalian preimplantation development.

Here, we present a modification of single-cell tagged reverse transcription protocol to study gene expression on a single-cell level or with limited RNA input. We describe different enzymes for reverse transcription and cDNA amplification, modified lysis buffer, and additional clean-up steps before cDNA amplification. We also detail an optimized single-cell RNA sequencing method for handpicked single cells, or tens to hundreds of cells, as input material to study mammalian preimplantation development. For complete details on the use and execution of this protocol, please refer to Ezer et al.1.

Long-chain polyphosphates inhibit type I interferon signaling and augment LPS-induced cytokine secretion in human leukocytes.

Inorganic polyphosphates are evolutionarily conserved bioactive phosphate polymers found as various chain lengths in all living organisms. In mammals, polyphosphates play a vital role in the regulation of cellular metabolism, coagulation, and inflammation. Long-chain polyphosphates are found along with endotoxins in pathogenic gram-negative bacteria and can participate in bacterial virulence. We aimed to investigate whether exogenously administered polyphosphates modulate human leukocyte function in vitro by treating the cells with 3 different chain lengths of polyphosphates (P14, P100, and P700). The long-chain polyphosphates, P700, had a remarkable capacity to downregulate type I interferon signaling dose dependently in THP1-Dual cells while only a slight elevation could be observed in the NF-κB pathway with the highest dose of P700. P700 treatment decreased lipopolysaccharide-induced IFNß transcription and secretion, reduced STAT1 phosphorylation, and downregulated subsequent interferon-stimulated gene expression in primary human peripheral blood mononuclear cells. P700 also augmented lipopolysaccharide-induced secretion of IL-1α, IL-1ß, IL-4, IL-5, IL-10, and IFNγ. Furthermore, P700 has previously been reported to increase the phosphorylation of several intracellular signaling mediators, such as AKT, mTOR, ERK, p38, GSK3α/ß, HSP27, and JNK pathway components, which was supported by our findings. Taken together, these observations demonstrate the extensive modulatory effects P700 has on cytokine signaling and the inhibitory effects specifically targeted to type I interferon signaling in human leukocytes.

Comprehensive characterization of the embryonic factor LEUTX.

The paired-like homeobox transcription factor LEUTX is expressed in human preimplantation embryos between the 4- and 8-cell stages, and then silenced in somatic tissues. To characterize the function of LEUTX, we performed a multiomic characterization of LEUTX using two proteomics methods and three genome-wide sequencing approaches. Our results show that LEUTX stably interacts with the EP300 and CBP histone acetyltransferases through its 9 amino acid transactivation domain (9aaTAD), as mutation of this domain abolishes the interactions. LEUTX targets genomic cis-regulatory sequences that overlap with repetitive elements, and through these elements it is suggested to regulate the expression of its downstream genes. We find LEUTX to be a transcriptional activator, upregulating several genes linked to preimplantation development as well as 8-cell-like markers, such as DPPA3 and ZNF280A. Our results support a role for LEUTX in preimplantation development as an enhancer binding protein and as a potent transcriptional activator.

Computational approach to evaluate scRNA-seq data quality and gene body coverage with SkewC.

SkewC is a single-cell RNA sequencing (scRNA-seq) data quality evaluation tool. The approach is based on determining gene body coverage, and its skewness, as a quality metric for each individual cell. SkewC distinguishes between two types of single cells: typical cells with prototypical gene body coverage profiles and skewed cells with skewed gene body coverage profiles. SkewC can be used on any scRNA-seq data as it is independent from the underlying technology used to generate the data. For complete details on the use and execution of this protocol, please refer to Abugessaisa et al. (2022).1.

B1 SINE-binding ZFP266 impedes mouse iPSC generation through suppression of chromatin opening mediated by reprogramming factors.

Induced pluripotent stem cell (iPSC) reprogramming is inefficient and understanding the molecular mechanisms underlying this inefficiency holds the key to successfully control cellular identity. Here, we report 24 reprogramming roadblock genes identified by CRISPR/Cas9-mediated genome-wide knockout (KO) screening. Of these, depletion of the predicted KRAB zinc finger protein (KRAB-ZFP) Zfp266 strongly and consistently enhances murine iPSC generation in several reprogramming settings, emerging as the most robust roadblock. We show that ZFP266 binds Short Interspersed Nuclear Elements (SINEs) adjacent to binding sites of pioneering factors, OCT4 (POU5F1), SOX2, and KLF4, and impedes chromatin opening. Replacing the KRAB co-suppressor with co-activator domains converts ZFP266 from an inhibitor to a potent facilitator of iPSC reprogramming. We propose that the SINE-KRAB-ZFP interaction is a critical regulator of chromatin accessibility at regulatory elements required for efficient cellular identity changes. In addition, this work serves as a resource to further illuminate molecular mechanisms hindering reprogramming.

The Continued Absence of Functional Germline Stem Cells in Adult Ovaries.

Ovaries are central to development, fertility, and reproduction of women. A particularly interesting feature of ovaries is their accelerated aging compared to other tissues, leading to loss of function far before other organs senesce. The limited pool of ovarian follicles is generated before birth and once exhausted, menopause will inevitably commence around the age of 50 years marking the end of fertility. Yet, there are reports suggesting the presence of germline stem cells and neo-oogenesis in adult human ovaries. These observations have fueled a long debate, created experimental fertility treatments, and opened business opportunities. Our recent analysis of cell types in the ovarian cortex of women of fertile age could not find evidence of germline stem cells. Like before, our work has been met with critique suggesting methodological shortcomings. We agree that excellence starts with methods and welcome discussion on the pros and cons of different protocols. In this commentary, we discuss the recent re-interpretation of our work.

ATMP-classified, scalable, autologous cell spray for the treatment of skin wounds and assessment of its effects on wound healing clinically and on a molecular level.

BACKGROUND: Autologous split-thickness skin grafts (STSGs) are the standard of care for closure of deep and large burns. However, perforation and extensive fishnet-like expansion of the grafts to achieve greater area wound coverage can lead to treatment failures or esthetically poor healing outcomes and scarring. The purpose of this study was to validate an autologous advanced therapy medicinal product (ATMP)-compliant skin cell suspension and evaluate its efficacy to promote epithelialization. METHODS: Cells isolated from a piece of STSG according to ATMP classification requirements were sprayed onto 20 patients during a single operation in a validation study. Comparative evaluation of treatment efficacy was carried out using side-by-side skin graft donor site wounds that were standardized in depth. Firstly, we characterized wound healing transcriptomes at 14 and 21 days from serial wound biopsies in seven patients. Then, side-by-side wounds in four patients were treated with or without the skin cells. The wounds were photographed, clinical outcomes assessed, and the treatment and control wound transcriptomes at 14 days were compared to the untreated wounds' healing transcriptomes. RESULTS: The average cell yield after isolation from the STSG was 2.4 × 106 cells/cm2 with 96 % viability. The product contained mainly keratinocytes and their precursors but also other skin cells such as fibroblasts were present. As compared to vehicle-treated donor site wounds, the wounds treated with cells demonstrated improved epithelialization by both direct comparison and machine learning analysis of the transcriptomes. CONCLUSIONS: We showed that rapid and scalable ATMP-classified processing of skin cells is feasible, and application of the skin cells effectively promotes healing and epithelization of donor site wounds.

Human endometrial cell-type-specific RNA sequencing provides new insights into the embryo-endometrium interplay.

STUDY QUESTION: Which genes regulate receptivity in the epithelial and stromal cellular compartments of the human endometrium, and which molecules are interacting in the implantation process between the blastocyst and the endometrial cells? SUMMARY ANSWER: A set of receptivity-specific genes in the endometrial epithelial and stromal cells was identified, and the role of galectins (LGALS1 and LGALS3), integrin ß1 (ITGB1), basigin (BSG) and osteopontin (SPP1) in embryo-endometrium dialogue among many other protein-protein interactions were highlighted. WHAT IS KNOWN ALREADY: The molecular dialogue taking place between the human embryo and the endometrium is poorly understood due to ethical and technical reasons, leaving human embryo implantation mostly uncharted. STUDY DESIGN SIZE DURATION: Paired pre-receptive and receptive phase endometrial tissue samples from 16 healthy women were used for RNA sequencing. Trophectoderm RNA sequences were from blastocysts. PARTICIPANTS/MATERIALS SETTING METHODS: Cell-type-specific RNA-seq analysis of freshly isolated endometrial epithelial and stromal cells using fluorescence-activated cell sorting (FACS) from 16 paired pre-receptive and receptive tissue samples was performed. Endometrial transcriptome data were further combined in silico with trophectodermal gene expression data from 466 single cells originating from 17 blastocysts to characterize the first steps of embryo implantation. We constructed a protein-protein interaction network between endometrial epithelial and embryonal trophectodermal cells, and between endometrial stromal and trophectodermal cells, thereby focusing on the very first phases of embryo implantation, and highlighting the molecules likely to be involved in the embryo apposition, attachment and invasion. MAIN RESULTS AND THE ROLE OF CHANCE: In total, 499 epithelial and 581 stromal genes were up-regulated in the receptive phase endometria when compared to pre-receptive samples. The constructed protein-protein interactions identified a complex network of 558 prioritized protein-protein interactions between trophectodermal, epithelial and stromal cells, which were grouped into clusters based on the function of the involved molecules. The role of galectins (LGALS1 and LGALS3), integrin ß1 (ITGB1), basigin (BSG) and osteopontin (SPP1) in the embryo implantation process were highlighted. LARGE SCALE DATA: RNA-seq data are available at www.ncbi.nlm.nih.gov/geo under accession number GSE97929. LIMITATIONS REASONS FOR CAUTION: Providing a static snap-shot of a dynamic process and the nature of prediction analysis is limited to the known interactions available in databases. Furthermore, the cell sorting technique used separated enriched epithelial cells and stromal cells but did not separate luminal from glandular epithelium. Also, the use of biopsies taken from non-pregnant women and using spare IVF embryos (due to ethical considerations) might miss some of the critical interactions characteristic of natural conception only. WIDER IMPLICATIONS OF THE FINDINGS: The findings of our study provide new insights into the molecular embryo-endometrium interplay in the first steps of implantation process in humans. Knowledge about the endometrial cell-type-specific molecules that coordinate successful implantation is vital for understanding human reproduction and the underlying causes of implantation failure and infertility. Our study results provide a useful resource for future reproductive research, allowing the exploration of unknown mechanisms of implantation. We envision that those studies will help to improve the understanding of the complex embryo implantation process, and hopefully generate new prognostic and diagnostic biomarkers and therapeutic approaches to target both infertility and fertility, in the form of new contraceptives. STUDY FUNDING/COMPETING INTERESTS: This research was funded by the Estonian Research Council (grant PRG1076); Horizon 2020 innovation grant (ERIN, grant no. EU952516); Enterprise Estonia (grant EU48695); the EU-FP7 Marie Curie Industry-Academia Partnerships and Pathways (IAPP, grant SARM, EU324509); Spanish Ministry of Economy, Industry and Competitiveness (MINECO) and European Regional Development Fund (FEDER) (grants RYC-2016-21199, ENDORE SAF2017-87526-R, and Endo-Map PID2021-127280OB-100); Programa Operativo FEDER Andalucía (B-CTS-500-UGR18; A-CTS-614-UGR20), Junta de Andalucía (PAIDI P20_00158); Margarita Salas program for the Requalification of the Spanish University system (UJAR01MS); the Knut and Alice Wallenberg Foundation (KAW 2015.0096); Swedish Research Council (2012-2844); and Sigrid Jusélius Foundation; Academy of Finland. A.S.-L. is funded by the Spanish Ministry of Science, Innovation and Universities (PRE2018-085440). K.G.-D. has received consulting fees and/or honoraria from RemovAid AS, Norway Bayer, MSD, Gedeon Richter, Mithra, Exeltis, MedinCell, Natural cycles, Exelgyn, Vifor, Organon, Campus Pharma and HRA-Pharma and NIH support to the institution; D.B. is an employee of IGENOMIX. The rest of the authors declare no conflict of interest.

Transient DUX4 expression in human embryonic stem cells induces blastomere-like expression program that is marked by SLC34A2.

Embryonic genome activation (EGA) is critical for embryonic development. However, our understanding of the regulatory mechanisms of human EGA is still incomplete. Human embryonic stem cells (hESCs) are an established model for studying developmental processes, but they resemble epiblast and are sub-optimal for modeling EGA. DUX4 regulates human EGA by inducing cleavage-stage-specific genes, while it also induces cell death. We report here that a short-pulsed expression of DUX4 in primed hESCs activates an EGA-like gene expression program in up to 17% of the cells, retaining cell viability. These DUX4-induced cells resembled eight-cell stage blastomeres and were named induced blastomere-like (iBM) cells. The iBM cells showed marked reduction of POU5F1 protein, as previously observed in mouse two-cell-like cells. Finally, the iBM cells were successfully enriched using an antibody against NaPi2b (SLC34A2), which is expressed in human blastomeres. The iBM cells provide an improved model system to study human EGA transcriptome.

DUX4 is a multifunctional factor priming human embryonic genome activation.

Double homeobox 4 (DUX4) is expressed at the early pre-implantation stage in human embryos. Here we show that induced human DUX4 expression substantially alters the chromatin accessibility of non-coding DNA and activates thousands of newly identified transcribed enhancer-like regions, preferentially located within ERVL-MaLR repeat elements. CRISPR activation of transcribed enhancers by C-terminal DUX4 motifs results in the increased expression of target embryonic genome activation (EGA) genes ZSCAN4 and KHDC1P1. We show that DUX4 is markedly enriched in human zygotes, followed by intense nuclear DUX4 localization preceding and coinciding with minor EGA. DUX4 knockdown in human zygotes led to changes in the EGA transcriptome but did not terminate the embryos. We also show that the DUX4 protein interacts with the Mediator complex via the C-terminal KIX binding motif. Our findings contribute to the understanding of DUX4 as a regulator of the non-coding genome.

Viral infection-related gene upregulation in monocytes in children with signs of ß-cell autoimmunity.

OBJECTIVE: The pathogenesis of type 1 diabetes (T1D) is associated with genetic predisposition and immunological changes during presymptomatic disease. Differences in immune cell subset numbers and phenotypes between T1D patients and healthy controls have been described; however, the role and function of these changes in the pathogenesis is still unclear. Here we aimed to analyze the transcriptomic landscapes of peripheral blood mononuclear cells (PBMCs) during presymptomatic disease. METHODS: Transcriptomic differences in PBMCs were compared between cases positive for islet autoantibodies and autoantibody negative controls (9 case-control pairs) and further in monocytes and lymphocytes separately in autoantibody positive subjects and control subjects (25 case-control pairs). RESULTS: No significant differential expression was found in either data set. However, when gene set enrichment analysis was performed, the gene sets "defence response to virus" (FDR <0.001, ranking 2), "response to virus" (FDR <0.001, ranking 3) and "response to type I interferon" (FDR = 0.002, ranking 12) were enriched in the upregulated genes among PBMCs in cases. Upon further analysis, this was also seen in monocytes in cases (FDR = 0.01, ranking 2; FDR = 0.04, ranking 3 and FDR = 0.02, ranking 1, respectively) but not in lymphocytes. CONCLUSION: Gene set enrichment analysis of children with T1D-associated autoimmunity revealed changes in pathways relevant for virus infection in PBMCs, particularly in monocytes. Virus infections have been repeatedly implicated in the pathogenesis of T1D. These results support the viral hypothesis by suggesting altered immune activation of viral immune pathways in monocytes during diabetes.

SkewC: Identifying cells with skewed gene body coverage in single-cell RNA sequencing data.

The analysis and interpretation of single-cell RNA sequencing (scRNA-seq) experiments are compromised by the presence of poor-quality cells. For meaningful analyses, such poor-quality cells should be excluded as they introduce noise in the data. We introduce SkewC, a quality-assessment tool, to identify skewed cells in scRNA-seq experiments. The tool's methodology is based on the assessment of gene coverage for each cell, and its skewness as a quality measure; the gene body coverage is a unique characteristic for each protocol, and different protocols yield highly different coverage profiles. This tool is designed to avoid misclustering or false clusters by identifying, isolating, and removing cells with skewed gene body coverage profiles. SkewC is capable of processing any type of scRNA-seq dataset, regardless of the protocol. We envision SkewC as a distinctive QC method to be incorporated into scRNA-seq QC processing to preclude the possibility of scRNA-seq data misinterpretation.

CRISPR activation enables high-fidelity reprogramming into human pluripotent stem cells.

Conventional reprogramming methods rely on the ectopic expression of transcription factors to reprogram somatic cells into induced pluripotent stem cells (iPSCs). The forced expression of transcription factors may lead to off-target gene activation and heterogeneous reprogramming, resulting in the emergence of alternative cell types and aberrant iPSCs. Activation of endogenous pluripotency factors by CRISPR activation (CRISPRa) can reduce this heterogeneity. Here, we describe a high-efficiency reprogramming of human somatic cells into iPSCs using optimized CRISPRa. Efficient reprogramming was dependent on the additional targeting of the embryo genome activation-enriched Alu-motif and the miR-302/367 locus. Single-cell transcriptome analysis revealed that the optimized CRISPRa reprogrammed cells more directly and specifically into the pluripotent state when compared to the conventional reprogramming method. These findings support the use of CRISPRa for high-quality pluripotent reprogramming of human cells.

Generation of RNA sequencing libraries for transcriptome analysis of globin-rich tissues of the domestic dog.

We have developed a protocol for barcoded cDNA libraries of 48 samples to study gene expression across tissues in the domestic dog, Canis familiaris, by modifying the Single-Cell Tagged Reverse Transcription (STRT) protocol (Islam et al., 2012, 2014). The cDNA reads represent mRNA 5' ends, enabling the study of transcription start sites (TSS). Our modifications include longer UMIs for molecular counting and Globin-Lock® to deplete globin mRNAs that are abundant in blood and blood-rich tissues dominating all reads.

Embryonic LTR retrotransposons supply promoter modules to somatic tissues.

Long terminal repeat (LTR) retrotransposons are widely distributed across the human genome. They have accumulated through retroviral integration into germline DNA and are latent genetic modules. Active LTR promoters are observed in germline cells; however, little is known about the mechanisms underlying their active transcription in somatic tissues. Here, by integrating our previous transcriptome data set with publicly available data sets, we show that the LTR families MLT2A1 and MLT2A2 are primarily expressed in human four-cell and eight-cell embryos and are also activated in some adult somatic tissues, particularly pineal gland. Three MLT2A elements function as the promoters and first exons of the protein-coding genes ABCE1, COL5A1, and GALNT13 specifically in the pineal gland of humans but not in that of macaques, suggesting that the exaptation of these LTRs as promoters occurred during recent primate evolution. This analysis provides insight into the possible transition from germline insertion to somatic expression of LTR retrotransposons.

Epitranscriptomics of Ischemic Heart Disease-The IHD-EPITRAN Study Design and Objectives.

Epitranscriptomic modifications in RNA can dramatically alter the way our genetic code is deciphered. Cells utilize these modifications not only to maintain physiological processes, but also to respond to extracellular cues and various stressors. Most often, adenosine residues in RNA are targeted, and result in modifications including methylation and deamination. Such modified residues as N-6-methyl-adenosine (m6A) and inosine, respectively, have been associated with cardiovascular diseases, and contribute to disease pathologies. The Ischemic Heart Disease Epitranscriptomics and Biomarkers (IHD-EPITRAN) study aims to provide a more comprehensive understanding to their nature and role in cardiovascular pathology. The study hypothesis is that pathological features of IHD are mirrored in the blood epitranscriptome. The IHD-EPITRAN study focuses on m6A and A-to-I modifications of RNA. Patients are recruited from four cohorts: (I) patients with IHD and myocardial infarction undergoing urgent revascularization; (II) patients with stable IHD undergoing coronary artery bypass grafting; (III) controls without coronary obstructions undergoing valve replacement due to aortic stenosis and (IV) controls with healthy coronaries verified by computed tomography. The abundance and distribution of m6A and A-to-I modifications in blood RNA are charted by quantitative and qualitative methods. Selected other modified nucleosides as well as IHD candidate protein and metabolic biomarkers are measured for reference. The results of the IHD-EPITRAN study can be expected to enable identification of epitranscriptomic IHD biomarker candidates and potential drug targets.